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Abstract 

We consider the possibility of an alternative gravity theory explaining the 
dynamics of galactic systems without dark matter. From very general as- 
sumptions about the structure of a relativistic gravity theory we derive a 
general expression for the metric to order (w/c)^. This allows us to compare 
the predictions of the theory with various experimental data: the Newtonian 
limit, light deflection and retardation, rotation of galaxies and gravitational 
lensing. Our general conclusion is that the possibility for any gravity theory 

to explain the behaviour of galaxies without dark matter is rather improbable. 
PACS numbers: 04.50.+h, 04.25.Nx, 95.35.+d 
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Einstein's gravity theory and some other alternative gravity models are in good agree- 
ment with the experimental data in the solar system and the laboratory [|T|. However, the 
behavior of galactic systems poses a great challenge to gravity theories. For virtually all 
spiral galaxies the tangential rotational velocity curves tend toward some constant value. 
This fact is in sharp contradiction with the visible star (luminosity) distribution and the 
laws of Newtonian dynamics. The stars in the outer parts of galaxies rotate several times 
faster than predicted by the standard gravity theory. A similar problem is observed in grav- 
itational lensing f^. Just as in the solar system problems of the past century (concerning 
the orbits of Uranus and Mercury), there are two ways to resolve these difficulties. 

One, the most widely adopted, is the dark matter hypothesis 0. It is presumed that the 
visible stars are imbedded in a massive nearly spherical halo of nonluminous matter. The 
mass of the halo varies from one galaxy to another but generally it constitutes about 90% of 
the total mass. This hypothesis explains the flat rotational curves of galaxies. Yet it has its 
own troubles, in particular (i) no good model for the formation of the dark halo is known, 
and (ii) after much effort and many proposals no known form of matter has yet given a 
satisfactory model for the massive halo. (The few recently observed cases of gravitational 
micro lensing |Q are as yet far from conclusive evidence for the dark matter explanation). 

The second way is to assume that for galactic distances Newton's gravity law is no longer 
valid. This possibility has also been the subject of some discussions In particular it 

was shown that a modified gravitational potential of the form 

V^ = ^^[l + «e-^/-], (1) 
r(l + a) L -I 

where a = — 0.9,ro ~ 30kpc can explain fiat rotational curves for most of the galaxies. 
The potential (|I]) differs from the usual one by an extra exponential term. For the solar 
system this term equals 1 with high accuracy and (|I|) reduces to the standard form ip = 
— GM/r but for distances significantly greater than 30kpc the exponential term vanishes 
and we have once again a Newtonian potential ip = —GM/{1 + a)r but now with an 
approximately ten times bigger gravitational constant G/{1 + a). Besides (P several other 



modified gravitational models were considered in the literature 0. Most were introduced 
purely phenomenologically without derivation from some gravity theory. All of these models 
and also various attempts to construct nonrelativistic gravity theories have the same 
trouble. They cannot describe the motion of light without additional assumptions. A 
description of light requires a relativistic gravity model. 

We investigate the general possibility of constructing a relativistic gravity theory which 
can explain galactic mysteries and other experimental data like the classical solar system 
tests [111. We formulate some very general postulates about the structure of the theory: 

(i) Gravitational phenomena are described by the metric of space-time g^v and possibly 
some other set of fields ^ a- The theory is invariant under general coordinate transforma- 
tions. 

(ii) The trajectories of (structureless) massive test particles and light are timelike and 
null geodesies of the metric g^u respectively. 

(iii) The source for the gravitational fields g^^ and ^ a are the energy-momentum tensor 
T^j^ and some current J a- For the solar system and galaxies these sources can be taken in 
the form T^^ = Tqq = p, = in the {v/cY approximation. 

(iv) The theory has a good linear approximation. 

(v) Flat spacetime gf^i, = rj^u = diag(— 1, 1, 1, 1) and "^a = (^a some constant or 
almost constant field) can be considered as the background field configuration for the solar 
system and galactic scales. 

(vi) The theory does not possess any unusual gauge freedom for the metric field besides 
general coordinate invariance. Any gauge freedom for the field \E'a is fixed by an appropriate 
gauge fixing condition. 

(vii) The theory is not a higher- derivative theory. 

Let us briefly discuss postulates (i)-(vii). The postulates (i),(ii) are the usual postulates 
of the so called metric theory of gravity [|T|. Postulate (iii), especially the condition Ja = 0, is 
of crucial importance for our study. First, it allows us to make definite predictions about the 
post-Newtonian approximation of the theory without needing detailed information about 
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its structure. Second, as we'll see below, it ensures that the theory does not violate the 
equivalence principle. The nature of the current J a may be different; for some models it 
may be absent explicitly. In particular, for the Brans-Dicke theory ^ a is the scalar 
field and it does not have any corresponding matter source. However, for the Poincare 
gauge theory of gravitation [|T^ is the space-time torsion and the current is the spin- 
tensor of matter which vanishes to a high approximation since both the solar system and 
the galaxies do not contain large amounts of spin-polarized matter. Postulate (iv) excludes 
from our considerations all essentially nonlinear theories for which the linear approximation 
is invalid. Postulate (v) means that we neglect global cosmological effects. Postulate (vi) 
ensures that gravitational equations are nondegenerate. The assumption about the absence 
of extra gauge freedom for the metric is quite natural since such invariance normally imposes 
unphysical constraints on the energy-momentum tensor of matter fields. Postulate (vii) 
excludes from our consideration theories with Green functions of the form 1 / (□— m^)", n > 2. 
These postulates are actually not too restrictive. For example they allow a large class 
of geometric gravity theories derived from Lagrangians depending on the metric and the 
connection through the torsion, curvature and non-metricity. 

We are going to compare the predictions of any gravity theory which satisfies postulates 
(i)-(vii) with the experimental data from the solar system and the galaxies. A key point 
is that all of these systems are essentially post-Newtonian slow-motion, weak-gravitational- 
field systems Hence we can consider our theory in the linearized approximation and we 
can use the small post-Newtonian parameter ip v"^ (we assume c = ^ = 1) in order to 
solve the gravitational equations approximately is a typical gravitational potential and v 
a typical velocity in the system). Observe, that GToo = Gp = 0{v^) and, therefore, the 
leading corrections for the g^i,, "^a ai'e 0(f ^). Thus, we can represent our fields in the form 

9f,u = Vt^u + V' = + V-A, (2) 

where h^^^ipA = 0{v'^). We are interested in computing only the metric since the test 
bodies are not sensitive to the other fields. It is well known ^ that the first post-Newtonian 
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correction for the equations of motion of massive test particles depends on /iqo only. This is 
just the Newtonian approximation and — |/ioo is the gravitational potential. But the post- 
Newtonian equations for light include both the /iqo and /ij^, i,k = 1,2,3 components of the 
metric (the gravitational potential alone cannot describe the motion of light!). 

The h^iy, ipA can be obtained from the linearized equations of the theory. The invertible 
linear operator V of these equations is constructed with the help of 77^,^, and At 
this stage we work with the weak field relativistic approximation. In the end we obtain 
the required post-Newtonian approximation by dropping all terms with time derivatives Oq 
and replacing □ by A. We assume that T) does not contain terms of the form \E'|^(9°. Such 
terms usually violate spatial isotropy, verified experimentally with rather high accuracy [|l| . 
Now we can use the spin projection operators |]lT| and decompose our equations on the 
independent spin sectors. In general the field h^j_y can contain contributions of 4 different 
kinds of particles of spin 2"*", 1~, 0"*", 0"*". The corresponding projectors read 



,2+ 



0+ 
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here u^^, = d^d^/D, and 9^^, = rj^^ — u^^. The subscripts 1,2 in the spin 0"^ sector label two 
different kinds of particles of this type. A similar spin decomposition exists for the field ipA 
but we are not interested in its detailed contents since the corresponding source Ja vanishes. 



Nonzero contribution to the metric for the source Tqo = p can come only from P^^T^,^ and 
Pi^Tfj,^. All other projectors produce terms which either vanish in the given approximation 
or can be eliminated by an appropriate choice of coordinate system. In general each spin 
sector can contain several different particles. The form of the linearized equations in the 
spin sector is 

/ 

MooPo MoiPoi ■■• 

MioPio MnPi ■■■ Pi^A = , (3) 



V 



/ V 



/ V 



where the coefficients Mik are scalar polynomial functions of the operator □ = dad" (here 
for simphcity we omit the superscript 2+ and denote P^* = Pq? the subscripts 1,2,. . .label 
different 2+ modes). The determinant of the matrix Mij has the form ]\{^ — m^Y. The 
constants nip play the role of "masses" for the propagating linearized modes of our alternative 
gravity theory. Due to the orthogonality and completeness properties of the spin projectors 
the solution of the equations @ can be obtained merely by calculating the inverse matrix 
Nik = We are interested to know Nqq only; its general form is Nqq = J2 <^p/ — f^p)- 

Finally we have to replace all operators !/(□ — m^) by 1/(A — m^) which leads to Yukawa 
exponential potentials. The analysis of the spin 0"*" sector is completely analogous. Hence, 
we obtain a general form for the metric 



^00 = -1 + 2 

gik = Sik\i + 2 



712 n2+no 

(o-Q + To)U + apUp + '^'Pq 

p=l g=n2+l 



_ "2 _ n2+no 

^ p=i ^ 



(J=n2+1 



(4) 



Here the constants (Jp.,Tp depend on the parameters of the concrete model and/or the con- 
stant field The constants ap and Tp represent the contributions of the spin 2^*" and spin 
0"*" modes respectively. The Newtonian potential U and exponential potentials Up 



U = G 



P' 



\x — x'\ 



,d^x\ Ur. = G I ^ ^d^x', 



\X — X' 



(5) 



correspond to the massless and massive modes with mass mp. Now all information about a 
particular gravity theory is packed into several constants Ti^akifnp. 

We can compare the metric with the standard metric of the parametrized post- 
Newtonian (PPN) formalism in the same O(f^) approximation 



5^00 



-1 + 2U, 



gik = Sik (1 + . 



(6) 



Here the coefficient of U in goo is fixed by the Newtonian limit and the experimental value 
for 7 = 1 ± 10^^ comes from light defiection and retardation experiments in the solar system 
The only essential difference between (^) and (P) is that the standard PPN formalism 



6 



does not take into account a possible contribution from massive modes. Of course, the 
influence of massive exponential potentials on the predictions of the gravity model depends 
on the concrete values of the masses rrip. If the mass of the mode is large enough, i.e., such 
that l/rrip is significantly less than 1cm, then the contribution of Up can not be observed in 
gravity experiments and our metric reduces effectively to (|^). On the other hand if 1/mp 
is larger than the typical size of a galaxy then Up ^ U for galactic and shorter distances 
and we are left once again with the metric (P). In principle l/m^ could be about the size 
of the Earth or the solar system but in this case experimental data impose very strong 
restrictions on the magnitude of the constants dp, [|T2|. For example if 1/mp ^ lO^^cm 



then cTp, Tp < 10^*^. 

Therefore, we have hopes to explain the dynamics of galaxies if 1/mp has an intermediate 
value significantly larger than the size of the solar system but not greater than the typical size 
of a galaxy (compare with (P). For our purposes it is sufficient to consider the simple case 
with two massive particles of spins 2+ and 0"*" with approximately equal masses mi ~ m2 ~ 
lO^^^eV. Thus, we have with four unknown constants ctq, Tq, cti, r2. The experimental 
data imposes constraints on these parameters. In the solar system Ui ^ U2 ^ U and we 
have 

(To + To + (Ti + r2 = 1, |(o-o + ai) - (tq + r2) = 1, (7) 

where the first condition ensures the correct Newtonian limit while the second follows from 
experiments with light. Now let us consider distances larger than 30kpc. For these distances 
Ui ^ U2 ^ 0. As was mentioned above in order to reproduce the fiat rotational curves of 
galaxies without dark matter as far as is known we only need the gravitational interaction 
to be approximately 10 times stronger for large distances 

(To + To ^ 10, i(To - To ^ 10. (8) 

Here the second conditions follow from gravitational lensing, since it is known that the 
observed lensing is in conformity with the predictions of Einstein's gravity with dark matter. 



Thus, we have once again an approximately ten times stronger effective couphng constant 
for hght. Solving one has 

40 10 

(^0 ~ y, ~ -y , (^1 ~ -12, rs ^ 3 (9) 

Any gravity model with the parameters should be in good correspondence with exper- 
imental data in the solar system and should explain the behaviour of stars and light in 
galaxies with reasonable accuracy. Probably if we include more spin 0"*" and 2+ particles of 
various masses we could fit more details of the galactic rotation curves This hypothet- 
ical model is not so simple: the gravity theory should include, besides the usual graviton, 
at least two extra massive very light particles [|13|. The most serious trouble with such a 
theory comes from the negative sign of tq and (Xi. In accordance with a general theorem 
||l^ the sign of these constants must be positive for even spin fields and negative for odd 
spin fields. A wrong sign results in a propagating mode carrying negative energy; this is 
considered unacceptable in a theory. There is still a small possibility of escaping this prob- 
lem. It has been suggested that if the theory contains several particles with the same spin 
then the conditions on the parameters might be weakened In our opinion the possibility 
of successfully matching the galactic rotation curves while avoiding negative energy modes 
seems remote. 

On the other hand, a negative sign for the coupling constant is natural for odd spin 
particles; initially it was suggested that the exponential term in (|l]) is mediated by a spin 
1 vector particle In our scheme one can reproduce such a contribution if and only if 
Ja 7^ 0. In particular, it was suggested that Ja may be proportional to baryonic charge [jl5l . 



However, the baryonic charge to mass ratio varies from one body to another. Therefore, 



such an interaction is no longer universal and violates the weak equivalence principle |T6 
which has been verified experimentally to a very high accuracy |1|]. Now the importance 
and fundamental nature of postulate (iii) becomes clear. It ensures the universality of the 
gravitational interaction. All models which satisfy this condition should not have trouble 
with the equivalence principle at least on the modern experimental level. 

8 



We have compared our model only with part of the available experimental data. Even 
if one overlooks the problem with the wrong sign of the coupling constants the model must 
describe correctly, in addition to the already considered effects, the perihelion shift of Mer- 
cury, the energy loss of the binary-pulsar, and cosmological observations. Comparison with 
these data cannot be performed with the help of linearized equations and requires a more 
detailed consideration for each particular gravity theory. Of course, such comparison prob- 
ably result in additional perhaps severe restrictions on the theory under consideration (see 
e.g. [|l^] which discuss cosmological restrictions). 



Although our scheme covers a large class of the theories yet we can not exclude the 
possibility of constructing a model with the desired properties which violates one of our 
postulates. We want to mention briefly some possibilities which have been proposed. 

1. It has been shown |jl8| that a cosmological constant A ^ lO^^^cm^^ is able to ex- 
plain the fiat rotational curves of galaxies. This value of A is ten times bigger than the 
cosmologically acceptable limit. 

2. Consideration of quantum corrections to the Newtonian potential may result in some 
additional logarithmic long range terms The applicability of these results to galactic 



distances is not obvious. 

3. A quadratic in curvature lagrangian can also produce modified gravitational potentials 
with extra long range terms 0. This model has higher derivatives equations and requires a 
traceless energy- momentum tensor. 

4. It is possible to construct an essentially nonlinear theory with a non-quadratic kinetic 
term in the Lagrangian. This model explains the dynamics of galaxies due to the nonlinear 
nature of the equations in the regime of small accelerations [|16|,|[ . It was observed that such 
a theory may have troubles with faster-than-light waves. 

Although all these models explain the rotation of galaxies they have their own weak 
points, moreover they must be capable of predicting correctly other gravitational effects. 
Therefore we conclude that the possibility of explaining galactic mysteries with the help of 
a modified gravity theory looks quite improbable. 
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